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Abstract 

We summarize and extend the correlative definitions and principles of abstract operators, 
discuss the relation between abstract operators and pseudo-differential operators, add sev- 
eral new algorithms, furthermore, develop the theory of partial differential equations with 
abstract operators, and then systematically expound the basic methods. By combining ab- 
stract operators with the Laplace transform, we can easily derive the explicit solution of initial 
value problem of linear higher-order partial differential equations for n-dimensional space, and 
establish the general theory of linear higher-order partial differential equations. 
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1 Introduction 

In 1960s, the general theory of linear partial differential equations made important progress by 
using the generalized function and its Fourier transform, further, the theory of pseudo-differential 
operators. The significance of the latter theory lies in introducing the concept of symbol, but due 
to the lack of awareness of the analytic continuous fundamental theorem (See pQ), the concept 
of abstract operators remained unclear. Without the theories of abstract operators, the operator 
method of differential equations is just a technique, and the solution of partial differential equations 
is normally complex. In 1997, the author delivered the analytic continuous fundamental theorem, 
based on which the author introduced the concept of abstract operators, and derived the algorithms 
of five types of abstract operators (See pQ) such as 

exp(hd x ), s'm(hd x ), cos(hd x ), s'mh(hd x ) and cosh(hd x ). 

Where cos(ihd x ) — cosh(hd x ), sin(ihd x ) = i smh(hd x ) . Obviously, the abstract operator uses the 
same notation as infinite order differential operator, as it is the extension of the latter. However, 
within the frame of infinite order differential operators, complex operators are difficult to define, 
such as 

exp(iP(cy), cosCiP^) 1 / 2 ), p{ Q x)1/2 - , cosh^P^) 1 / 2 ), p {dx y/2 ■ 
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But as abstract operators, we can easily establish their algorithms (See [T], [2]), therefore the general 
solving procedure of initial value problem of linear partial differential equations with constant 
coefficients is derived pQ , [2] , [3] , this approach can also be applied in solving the partial differential 
equations with variable coefficients containing t. 

According to [4], pseudo-differential operators are defined as: 

If a(x,g) € C°°(IR™ x R£), and for arbitrary a,/3 G N n and an real number m, 

\d^a(x,o\<c a Ai + mr- lal 

is tenable, where C a< p is a constant, then the linear continuous mapping A of S"(K n ) — > ^(W 1 ) 
can be defined as: 

Au(x) = (2n)- n J e^a(x, OHtH, (I) 

which are called the pseudo-differential operators, denoted by a(x, D), where a(x, £ ) are the symbols 
of a(x, D), D denotes 

Obviously, this definition means the following equation is tenable, namely 

a(x,D)e^ x = a{x,Oe^ x . (2) 

Therefore pseudo-differential operators are similar to abstract operators. As the definition formula 
of pseudo-differential operators is a Fourier integral expression, their symbols are largely restricted 
to ensure the convergence of Fourier integrals. But for many specific problems, the symbols of 
pseudo-differential operators should not be restricted that much. Thus, a better definition of 
symbols can be expressed as: 

a A {x^) = e~ lix Ae* x . (3) 

Where A are properly extended pseudo-differential operators. Certainly, this extension still has its 
limitations as it involves another type of Fourier integrals. 

Then can we get rid of Fourier integrals when defining pseudo-differential operators? Let us 
study the following expression: 

Ae** =a A (x,Z)e iix . (4) 

Obviously, the properties of A depend on va{x,£,)- (0]) indicates the mapping relation between A 
and oa{x,£,)- For a given function 17,4(2;, £), if the corresponding algorithms of A can be derived 
from this type of mapping relations, to determine the domain and range of A, then this type of 
mapping relations can be the best definition of pseudo-differential operators. Luckily, this idea 
totally works by using the analytic continuous fundamental theorem or the abstract operator 
fundamental theorem. The pseudo-differential operators established by using this method are the 
abstract operators. 

Now let us carry out a systematic review of the basic methods of abstract operators, and prove 
that if we combine abstract operators with the Laplace transform on C°°(R™ xRj), then no matter 
how complex the initial value problem of linear partial differential equations is, the solving process 
is simply easy. In addition, in order to achieve self consistency theoretically, it is necessary to 
derive rules of differentiation from the view of abstract operators. 
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2 Basic theory of abstract operators 



2.1 Basic concepts, research method and basic operations 

Definition 1. (See pQ) Within each convergence circle of analytic functions, if the effects on 
the series term by term from the linear operator converge uniformly to the effects on the sum 
function, then the operator is called having the analytic continuity. 

Definition 2. The linear continuous mapping A of J^(R n ) — > J^(K") has the Fourier continu- 
ity, if A acts on an arbitrary function / £ ^(R n ), we have 

Where the Fourier transform / = & f is defined as 

The Fourier transform & is an isomorphic mapping in the space ^(W 1 ), and its inverse transform 
is 

Definition 3. (See pQ) a: Q , a; S R™ is called the base function, its exponential form is e^ x , 
where a £ N™, £ £ R n or C" are respectively called the character of x a and e* a . 

Definition 4. (See pQ) A class of linear operators is called the abstract operators, denoted by 
f(t, d x ), x £ K n , t £ R 1 , and if it acts on the exponential base function e^ x , we have 

f(t,d x )e& = f(t,Z)e&, V/(t,0 £ C°°(B} x R£). 

Which is called the abstract operators taking d x as the operator element, where f(t,£) is called 
the symbol of the abstract operators f(t,d x ). 

If abstract operators f(t, d x ) is a linear operator having the Fourier continuity, then it similar 
to a pseudo-differential operators. 

Definition 5. A class of linear operators is called the abstract operators, denoted by / (t, , 
x £ R", t £ R 1 , where 

8 _f d d d \ 

aa; \ ctei 0x2 ox n J 

and if it acts on the base function x a , we have 

/ U x-^j x a = f(t, a)x a , V/(t, a) G C 00 ^ 1 x N"). 

Which is called the abstract operators taking as the operator element, where f(t,a) is called 
the symbol of the abstract operators / (t,x-^). 

Definition 6. We call any operator element and common variables the abstract element, 
denoted by X,Y,... and X = (A 1; X 2 , • . . , X n ), Y = (Y ll Y 2 , . ■ . , Y n ), etc.; Accordingly, f(X), 
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f(Y), etc., represent the abstract operators or common functions, which is called the operator 
function. 

Abstract operators fundamental theorem. Let A,B, A',B' be the abstract operators, if 
there are functions v(t) £ 9(A'), u(t) € 9(B'), teR 1 and abstract elements X,Y, making one 
of the following two expressions tenable 

AX a A'v(t) = BY a B'u{t) or Ae^ x A'v(t) = Be^ Y B'u(t), 

and the expressions of A, B, A',B' do not explicitly contain the character a E N™, £ € K„orC n , 
then 

Af(X)A'v(t) = Bf(Y)B'u(t), V/(X), f(Y) e 9(A) n 9(B). 

Where 9(-), &(•) are domains and ranges of operators. 9(f(X)) = M(A'), 9(f(Y)) = M(B'), 
9(A)=M(f(X)), 9(B)=M(f(Y)). 

When A'v(t) = B'u(t) = I, and X = x e R n , Y = y e W 1 , this theorem becomes the following 
Analytic continuous fundamental theorem: 

If there are functions x(t),y(t) € K", t € K 1 , making one of the following two expressions 
tenable 

Ax a =By a or Ae^ x = Be iy , 

and the expressions of abstract operators A, B do not explicitly contain the character a € N™, 
£ G R„ orC„, then 

Af(x) = Bf(y), Vf(x)J(y) G n 9(B). 

Especially, when A, B are the linear operators having the analytic continuity, the analytic contin- 
uous fundamental theorem has been proved by G.Q Bi in 1997 (See [1]). 

The abstract operators fundamental theorem is the fundamental property of abstract operators, 
which directly extends the domains of abstract operators from power functions and exponential 
functions to a wider range of functions. 

Definition 7. A linear operator equation, if it's true for arbitrary functions or operator 
functions in a certain range, then it is called an operator formula. Particularly, the operator 
formulas that determine the domain and range of abstract operators are called the algorithm of 
the abstract operators. 

Definition 8. The relational expression between each component of the character a = 
(a\,cx2, ■ ■ • ,a n ) or £ = (£i, £2, ■ • • , £n) is called the characteristic equation. In general, an op- 
erator formula becomes the corresponding characteristic equation, when where arbitrary function 
of which is the base function. 

Corollary 1. As the base function can be expressed by both power function and exponential 
function, a single operator formula can be corresponding to two different characteristic equations, 
conversely, a single characteristic equation can also be corresponding to two different operator 
formulas. 

Constructing the operator equation tenable for the base function by the characteristic equation 
through the definition of abstract operators, and then deducing that it is also tenable for arbitrary 
functions or operator functions in a certain range by using the analytic continuous fundamental 
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theorem or the abstract operators fundamental theorem, thus we derive new operator formulas. 
Finding or establishing a new operator formula requires the knowledge of the corresponding char- 
acteristic equation in advance, without knowing the specific form of the new operator formula. 
Therefore, it all boils down to seek or construct appropriate characteristic equations. The key 
to transform the characteristic equation to the corresponding operator formula, is constructing 
the operator equation true for the base function by using the specific form of the characteristic 
equation and the definition of abstract operators, and also, only when the operator constructed 
is linear and the expression of the linear operator doesn't explicitly contain the character of the 
base function can we derive that the operator equation is not only true for the base function, but 
also for arbitrary functions in a certain range, according to the analytic continuous fundamental 
theorem or the abstract operators fundamental theorem. This is the research method of abstract 
operators. 

Example 1. Let x G R n , £ G C R„, then V/(£) G C°°(fi) we have 



f(O x ) S inbx = sin(6x+^ 1 A + ^ 2 ^ + ... + |^|_ )/(0 

7T d 7T d 7T 9 

f(d x )cosbx = cos ( bx + - £1— + - £ 2 — + ••• + - £„— ) /(£) 



(5) 



Where 6a; = b%xi + 62^2 + • • • + b n x n . 

Proof. Selecting the following two expressions from the known formulas as the characteristic 
equations, namely 



d" sin bx — b a sin (bx + — \a\j and d™ cos bx = b a cos (bx + — 



a 



Taking a G N n as the character of the base function, d x and as the operator elements, we have 
the characteristic equation with the following two operator equations by Definition 5: 



tfdnta = sin ( bx+ | fl A + + ... + ^ n JL ) C 

7T 9 7T d 7T d 

d°COSbx = COS ( bx+ J &o7- + j6o7- + --- + JCnflT- ) C 



C=6 



2^96 2~0£ a 2^ n/ 5=b 

According to the abstract operators fundamental theorem, we get the Example 1. 

Definition 9. Let A be a linear operator having the analytic continuity or Fourier continuity, 
if there is another linear operator, denoted by A^ 1 , making AA^ 1 = A~ x A = I, then A^ 1 is called 
the inverse operator of A. 

By Definition 4 we have the following Corollary: 

Corollary 2. (See pQ) The operator algebras formed by all the abstract operators f(t,d x ), 
are isomorphic to the algebras formed by their symbols /(£,£)• This isomorphism is determined 
by f(t,d x ) O f(t,0, and 

/(*, d x ) ± g(t, d x ) O f(t, ± g(t, 0, f(t, d x ) o g(t, 8 X ) o /(t, 0- 

Especially, the abstract operators /(9a) and g(d x ), x G R™ are each other's inverse operator, if 
and only if their symbols /(£) and <?(£) satisfy = 1, £ G 1„. 
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For instance, d Xi su\(xd y ) = d Vi cos(xd y ), d Xi cos(xd y ) = —d Vi sm{xd v ), where x,y E M". 
Corollary 3. (See pQ) Let g{x) € C°°(R n ), if g{d^)(e^ x f(0) is continuous at £ = a, then 

f(d x )(e a *g(x)) = g(d s )(e^f(0)k= a . (6) 

Obviously, by Definition 5 we can also get similar Corollary. 
Example 2. If function f(x),x G R" satisfies 

f{rj) exp ^- ^ ~ 4 ^ X * ^ drj < +oo, r? e E n , a > 0, 
then 

1 f ( \r) — 2ax 12 



Proof. In the Corollary 3, let g(x) = e a ^ 2 E C°°(R") and a = 0, we have 
/(a,)^*! 2 = exp(a|^| 2 )( e ^/(O)| ?=0 = exp(a|x + ^| 2 )/(O| ?=0 
= e Q l-l 2 exp(2axa c )exp(a|a ? | 2 )/(O| ?=0 

a|x|' 1 —/o.._ a \ /• 1^7 ~ CI 2 

2™(a7r)™/ 2 



exp(2ax^) /fa) exp LJ^JL^ <fy 



/0?) exp dr] 



2™(a7r)™/2 J Rn jy " * \ 4a 



a > 0. 



Where xd^ — x\d^ + X2d^ 2 + • • • + x n d^ n , \x\ = xf + x% + ■ 

Clearly, by ([7]) the /(r?) can also be a bounded measurable function. Therefore, Corollary 3 is 
an important way to extend the symbols and domains of abstract operators. In fact, limitations 
on the symbol f(t,£) and the operator element X will be broadened, and more general abstract 
operators will be introduced to deal with problems involved in our further discussion. 

There are two major research methods of the operator theory. One determines the domain and 
range of an operator first, then derives its operator algebra. The other abstractly determines an 
operator algebra, then derives the domain and range of the operator. The former is applied by 
pseudo-differential operators, while the latter is applied by abstract operators. The definition of 
abstract operators actually determines a kind of excellent operator algebras, then establishes dif- 
ferent algorithms for different operators using the abstract operator fundamental theorem, further, 
derives different domains and ranges for abstract operators with different symbols. In other words, 
the domain and range of an abstract operator depend on the algorithm it applies. Therefore, a 
significant feature of abstract operators is focusing on the research of algorithms, without depend- 
ing on the Fourier integral transform. The basic theory of abstract operators is a theory about 
algorithms, which has already become an important tool to seek and calculate explicit solutions of 
linear higher-order partial differential equations. 

2.2 Algorithms of Abstract operators on different iable function 

In order to achieve theoretical rigor, firstly let us derive the rules of differentiation from the 
view of abstract operators, making the concept of ordinary or partial differential operator on the 
basis of the abstract operators itself. 
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Differential operators can be defined as the following abstract operators, namely 

x-^-x a =ax a , ^-e^ x = ^ x , VieK 1 . 
ax ax 

Where a € N 1 , £ G Ri are the symbols of abstract operators x-^ and 4- respectively. 

Let a and 6 be the characters of base functions, as the characteristic equation, the binomial 
formula of integer power can be expressed as 



(a + b) n = a?b n -i , Vn G N 



By base functions e ax and e bx we can combine the binomial formula with the following operator 
equations: 

( e ax . e bx\ = f n \ ^L_ e ax . dn ~ J c bx 

dx n y ' dx? dx n -i 

According to the analytic continuous fundamental theorem, it is tenable, namely 

^M-Ef^?? 2 ^ W.ueC-CR 1 ). (8) 

This is the Leibniz rule. Similarly we have 

«? f =E(- 1 ) i f n> l^f^V Vv,ueC«(^). (9) 

dx n j-^ ' \j J dx n -3 \ dxi ) 
By (fTCTj) we easily obtain the Hermite polynomials 

The major parts of rules of differentiation are the derivative principle of function product, the 
derivative principle of compound function and the chain rule of multivariate function. Firstly let 
us establish the derivative principle of function product. Without losing the universality, we only 
consider the situation of the product of two functions. Let n = 1 in ijSJ), then for v = f(x) and 
u = g(x) we have 

±(f( X )g( X )) = f( X )±g( X )+g( X )±f( X ), V/, g € C 1 (M 1 ). (11) 

Let ip(x) = fi(x)f 2 {x) ■ ■ ■ f n (x), generally we have 

j n n j 

3=1 K 
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Let fi(x) — f2(x) = ■ ■ ■ = ,fn{x) = g(x) = y, then for differentiable function y = g{x) we have 

d n n-i d V d n dy d n , 

—y = ny -j- or —y = -j-^-y , Vn € N . 

dx dx dx dx dy 

According to the analytic continuous fundamental theorem, we have the derivative principle of 
compound function f(g(x)) € C 1 (R 1 ): 

By using (ITT1) and (fT2)l . for differentiable functions x\(t) and x% (t), t € R 1 we have 

dt [Xl X2 j " din dt X2 +Xl dx 2 dt' ai,a2£A 

Taking this one as the characteristic equation, by base functions x a 7 x € R 2 , a £ N 2 , combining 
it with the following operator equation: 

d a= dx 1 d a+ dx 1 d a ^ xWeR 2 feR l aeN 2^ 

dt dt ox\ dt 0x2 

According to the analytic continuous fundamental theorem, it is tenable V/(x) £ C 1 (R 2 ), namely 
d „, , dxi d „. . (ixo 9 „, , , , m o m i 

Clearly, for £R n ,t£ R 1 , we can generally derive V/(x) G C x (M n ) 

1 f( x ) = ^lA/W + /(. T ) + . . . + ^ -^-f( x ). (13) 

dt (it dx\ dt dx 2 dt dx n 

If taking n G N as the character of base function, the binomial formula can be expressed as the 
following characteristic equation: 

(x + h) n = y %^-x n . 

3=0 J 

According to the analytic continuous fundamental theorem, it is tenable for any analytic function 
/(i)er(C 1 ), namely 

hi di 

/(x + f{x), x£C\ (14) 

j=0 ^ x 

This is the Taylor formula. If taking (|14[) as a characteristic equation, then similarly we have 

Theorem 1. If 3v(s),u(s) £ C 00 ^ 1 ) can make the infinite series on the right side of (TT51) 
uniform convergent, then it will uniform converges to the left side of the formula, namely 

/(4)<™>=|£S^4)»> ■^■wf « i5 » 

Corollary 4. In (|T5j) . when A = —1, v — s, u — 1/s, s ^ 0, we have 

f '(-§?)l =sf (-lL)-s- m ° r ^/'(*) = ^/(*)-/(0)- (16) 
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Where JSf/(i) denotes 

^f(t)=f(--p\-, > 0, « G Hi. (17) 



ds J s 

This is the internal connection between abstract operators and Laplace transform, thus the abstract 
operator becomes a significant tool to compute Laplace transform. For instance, by (|17p we have 

Jt?[gW(t)]=g(-^)F(s), (F(s)=J?f(t)). (18) 



Similar to (fTT|) - ([T^|) . Guangqing Bi has already established a series of algorithms of abstract 
operators in reference pQ: 

cos(hd x )f(x)=^[f(x + ih)}, sin(hd x )f(x) = $[f(x + ih)}, Vf(z)eC°°(C n ). (19) 

Where bx = b\X\ + 62X2 + • • • + b n x n , bh = b\h\ + 62^2 + • ■ ■ + b n h n . 



sm(hd x )(uv) — cos(hd x )v ■ sin(hd x )u + sm(hd x )v ■ cos(hd x )u, 
cos(hd x )(uv) = cos(hd x )v ■ cos(hd x )u — sm.(hd x )v ■ sm(hd x )u. 



■ /'j.a \ u cos(hd x )v ■ sm(hd x )u — sm(hd x )v ■ cos(hd x )u 



v (cos(hd x )v) 2 + (sm(hd x )v) 



2 



cos(hd x )- 



u 



sin (/>£)/(.(*)) = dn(^)/W, 
cos (h£\ f(x(t)) = cos(Y^f(X). 

Where Y = sm(hd t )x(t), X = cos(hd t )x(t), t e R 1 , h £ Hi orC x . 
If n = 2, then (|22l) can easily be restated as 



^(h^)f(x(t),y(t)) = sin(Y x -JL + Y y ^)f(X x ,X y ), 



cos (h^)f(x(t),y(t)) = C0S ( Y -a^+ y ^)^'^)- 

Where = sm(hd t )x(t), X x — cos(hd t )x(t), Y y = sm(hd t )y(t), X y — cos(hd t )y(t). 
Theorem 7 in reference pQ can be more generally expressed as: 



(20) 



cos(hd x )v ■ cos(hd x )u + sm(hd x )v ■ sm(hd x )u 
v (cos(hd x )v) 2 + (sin(hd x )v) 2 

Where x € R", h E R n orC„. 

sm( hd t )f(x(t)) = sm(Ydx)f(X), 
cos(hd t )f(x(t)) = cos(rax)/(X). 

Where X = (X 1 ,X 2 , ■ ■ ■ ,X n ), X 3 = cos(/i9 t )x i (t), Y G R„, Y,- = sin(/i9 t )x i (t). 
For example, when n = 1, the (|22|) can easily be restated as 



(23) 



(24) 
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Theorem 2. If y = f(bx) in set of analytic functions is the inverse function of bx = g(y), 
namely g(f(bx)) = bx, then sm(hd x )f(bx) (denoted by Y) and cos(/i<9a;)/(&ir)(denoted by X) can 
be determined by the following set of equations: 

C0S ( Y 1^)9(X) = bx, 



dX 

d \ x £ K n , b £ R n , h £ R„. (25) 



sm\Y—jg(X) = bh. 

By using (l20|) - (l25t . we can constructively extend the domains and ranges of the following five 
abstract operators 

exp(hd x ), sin(hd x ), cos(hd x ), sinh(hd x ) and cosh(hd x ) 
from power functions and exponential functions to all the C°° functions. For instance (See [5]) 



s'm(hd x ) arctan&a; = — tanh l - 



2bh 



1 + (bx) 2 + (bh) 2 ' 

1 26x (26) 
cos(no x ) arctanoir. = - arctan — — . 

2 1 — (ox) 2 — (bhy 

Therefore, establishing such a set of algorithms as (j20|) - (l25l is the second important way to extend 
the domains and ranges of abstract operators. 

Because of the introduction of Definition 5, it is necessary to add new algorithms. 

Theorem 3. Let a £ C 1 , p £ K 1 , then V/(z) £ C 00 ^ 1 ), we have 

a p ^ f(p) = f(ap). (27) 
Proof. By Definition 5, for a £ C 1 , pGl 1 , we have 



p ^ n a — n a n a 



or °p p = a 



p a or a p ^p a = (ap) a , Va £ N 1 . 



According to the analytic continuous fundamental theorem, we have Theorem 3. 

Theorem 4. Let X = (pi cos &\ , ■ ■ ■ , p n cos n ), Y = (p\ sin 0%, ■ ■ ■ , p n sin 9 n ) , introducing the 
notation 

d d d d d 

dp^ r = lPl — + --- + 6 nPn — , Yd x = Yi-— + --- + Y n - 



Op Opi Op n OXi dX„ 

then V/(z) £ C°°(C"), p £ M™ , 9 £ K„, we have 

cos ( 9 P 4-)f(p) = cos(Yd x )f(X), 



dp 

sin (8p-^f(p) = sm(Yd x )f(X). 
Proof. Let cii = e 1 ® 1 , . . . , a n = e n , according to Theorem 3, then 

exp (iOp-^j f(p) = f(pe ie ) = f(X + iY) = cxp(iYd x ).f(X), X £ R", Y £ M„. 



(28) 
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d \ /„ d 



Considering 

exp\iVp—J = cos\Vp—j+ism\Vp—.. 

exp(iYdx) = cos(Ydx) + i sin(Ydx), 

thus we have Theorem 4, which transforms the abstract operators taking p-j^ as the operator 
element into those taking dx as the operator element. 

Corollary 5. Let f(z), z = x + iy € C 1 be an arbitrary analytic function in a complex plane, 
then on the line y = kx we have 



cos [y^J f(x) 
sin (y^j f(x) 



= /3 X 9* cos I ax^- ) f(x), 



y—kx 



y—kx 



fi x L s i n [ ax-^- ) f(x) 



dx t 

d \ (29) 



Where j3 = (1 + fc 2 ) 1 / 2 , a = arctanfc, which can be used to solve boundary value problems of 
2-dimensional Laplace equation on polygonal domains. 

Corollary 6. Making use of the analytic continuous fundamental theorem and the Pochham- 
mer symbol (A) rn defined by 

_ r(A + m) _ [ 1, (m = 0), 

{)m ~ T(A) |A(A + l)-(A + m-l) (mGN+). 1 J 

It is easily seen 

(-x) m t l (i) = (-i) n i ra ^ 1 v^en^U^r^r^^- ( 31 ) 

Where x = e z or z — In a;. Using this set of formulas, we can easily transform the Euler equation 
into an ordinary differential equation with constant coefficients. 

Therefore, for an arbitrary abstract operators, its symbol are C°° functions. Further descrip- 
tions require concrete analysis under specific circumstances. For instance, the symbol of the ab- 
stract operators on the right side of (|17l) can be further described by using conditions of the Laplace 
transform. The domain of the abstract operators P(d x ), sm(hd x ) and cos(hd x ) are the C°° whole 
space. 

2.3 Abstract operators and bounded function 

The abstract operators can also acting on the differentiable functions on bounded domain. 
For instance, let a < b, without losing the universality, assuming a > 0, b > we have 

na' 1 - 1 < (a"- 1 + a n ~ 2 b + a n - 3 b 2 + ■■■+ ab 71 ' 2 + b 71 ' 1 ) < nb n ^\ 

n i b n -a n ,„ , /16™-a n \^ , 

na™" 1 < — < nb 71 - 1 or a < ( ; ) < b. 

n b — a 
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Therefore, if a < b, then we have a < c < b, making 



ric 



n-l 



b n - a n 



b n - a n 



b — a ac b — a 

Taking this one as the characteristic equation, according to the analytic continuous fundamental 
theorem we get the Lagrange mean value theorem, namely 
If a < b, and f(x) £ C 1 [a, b], then 3c £ [a,b] making 



Vn £ N 1 . 



d f( \ f( b ) - /(«) w,, \ c n l, .1 

— /(c) = , Vf(x)£C[a,b\. 

ac b — a 

In reference [5], Guangqing Bi and Yuekai Bi discussed the following functions: 



(32) 



f(x) = cos(-#)5( e ~ 



g(x) = sin ( Z__i_ ) S(e z ) 



c dz 
rx d 
c dz 



z=0 



z=0 



z £ M 1 , a < x < b. 



(33) 



Where f(x),g(x) £ C^R 1 ), S(t) £ C°°(r2 ), the ft is a certain neighborhood at t = 0. 

If 5(t) is expanded in power series on Qo, then (I33|) are equivalent to Fourier series, and their 
domain (a,b) can be extended to all (— oo,+oo) periodically. So f(x),g(x) £ L 2 {[— c, c]) and 
/(x + 2/) = /(x), s (i + 20=fl(a;). 

For instance, by (j23j) and (|26)) . we have 

d 



cos 



7TX 9 

— I arctan e* 



c 9z 



= cos Y 



z=0 



dX 



arctan 



1 2X 

2 arctan i-(^ + r 2) 



z=0 
2 COs(7Tx/c) 



2 = 



1 

- arctan -, 

2 1 — (cos 2 (7ra;/c) + sin (irx/c)) 



When cos(7ra/c) ^ or |x| ^ fee + c/2, fc = 0, 1, 2, • • •, the above expression turns into 



cos 



d 



dz 



j arctan e z 


z=0 \ 







(1/2) arctan(+oo) = +tt/4, 2Kc - c/2 < x < 2Kc + c/2 
(1/2) arctan(-oo) = -tt/4, 2Kc + c/2 < x < 2Kc + 3c/2. 

Where K = 0, ±1, ±2, • • •. This function represents a square wave. 

By algorithms (fl9 ]) -(|25 ]l . and some basic formulas similar to ((26)) (See [5]), we obtained 



sin 



7TX d 

c dz 
kx d 
c dz 



(-ln(l-e 2 )) 



«»( Mi -e*)) 



7TX 5 

sin| 7^ )ln(1 + r) 

cos ( — «I )ln(l + e z ) 



c <9,z 
c <9z 



z=0 



z=0 



z=0 



arctan 



7T 7TX 

2 ~ 2c' 



— in 2 sm — 

2c 



< x < 2c. 

< x < 2c. 



■KX 

2c"' 



\x\ < c. 



, , TTX 

In 2 cos — 

2c 



7T 

4' 



\x\ < c/2. 



\x\ < c. 



(34) 



In other words, the ranges of the abstract operators sm(hd x ) and cos(hd x ) may be extended from 
the QfOR 1 ) to the L 2 ([-c,c\). 
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In order to solve an initial-boundary value problems of partial differential equations, we need to 
express a function with boundary conditions. For instance, (|34p are infinitely differentiable within 
a finite domain, such as 



d k 




7ra;\ 


z=0 dxk 


(f 


~ 2c) 



, VfceN, < x < 2c. 



Example 3. V/(z) e C 00 ^ 1 ), if /(z) is continuous at z = 1/2, then 



<!■>■ = J ( \ 



(35) 



p— cos a: 

Proof. Taking the following integral formula as the characteristic equation 



cos™ x cos na; da; 



2 n+l 



Vn € N. 



Combining the characteristic equation with the following operator equation by using Definition 5 
and the base function x a : 



-/2 / q 

cos I x P-Q- p 



p n dx 



VneN 1 , xg[0,tt/2]. 



Where p = cosx, according to the analytic continuous fundamental theorem, we get (|35[) . 
Example 4. V/(z) € C^fjC 1 ), if f(z) is differentiable at z = 1/2, then 



ii 



*/2 / g 

8111 ^ 



dx 



p— cos a; 

Proof. Considering another integral formula 

r7r/2 

cos™ x sin nx dx = - 



1 /• 1 /(Q-/(l/2) 
1/2 



d£. 



(36) 



1 n 



Vn G N. 



fe=i 



Obviously it cannot be taken as the characteristic equation directly, but considering J2k=i 2 k /k = 
Jo ^q-i dq, and according to the analytic continuous fundamental theorem, we get (|36p . 

Example 5. Let k E N + , V/(z) £ C 00 ^ 1 ), if /(z) is /c times continuously differentiable at 
z = 1/2, then 

*7r/2 



cos('2/,\r) cos ( xpj- ) /(/;) 



2k+i kl 



2 r /2 ( d\ 

— J sm(2kx) sin \x p—\ f(p) 



dx 



— s {k) (- 



p— cos sc 



2 k + 1 kl 



Proof. In [B], Connon has given the following useful identical equations(Vn g N): 

2™ cos™ xcos(nx) — J cos(2fcr). 



fe=0 



(37) 
(38) 

(39) 
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2™ cos™ x sin(nx) = ^ ( j sin(2fcx) 
fc=o ^ ' 



Using the orthogonality of trigonometric functions, by ((39)) and (|40| we have 



i r 

— I cos(2fcx) cos™ x cos(nx)dx 



sin(2fcx) cos™ x sin(nx)<ix 



7T 

i r 

I" 7o 



fc) (2 

n\ / 1 



n+l 



, Vfc g N+. 



fc/ \2 



n+l 



, Vfc e N+. 



Which can also be written as 

2 r /2 



cos(2fcx) cos™ x cos(nx)dx 



tt/2 



sin(2fcx) cos™ x sin(nx)<ix = 



:) G 

1) a 



n+l 



n+l 



(40) 

(41) 
(42) 

(43) 
(44) 



Clearly, (|4"3")l is the characteristic equation of (|37l) , and |4"4")) is the characteristic equation of 
So by Definition 5 we have 

2 r /2 



vt»l 2k. c) * 'os ( x P-^~ ) />'" 



2 f^ 2 / <9 \ 

— / sin(2fcx) sin I x p— ) p r ' 
Wo V d P J 



dx 



1 d fc 
2 fe +!fc! d^ 



, VneM 1 . (45) 



{=1/2 



p=cos a: 



1 d fc 

dx = - — r 



, VneN 1 . (46) 



5=1/2 



Then according to the analytic continuous fundamental theorem, we get (|37[) and 
In reference [3] we get 

Let S(t) = J2n=o a » r ' ^R 1 , < i < r, < r < +00, if 



E a ™ 



nirx 
cos = cos 

c 



7TX 9 

c 9^ 



E 

n=0 



totx / nx d 

a n sin = sin — — 

c \ c az 

b-o 



S(e z ) 
S(e z ) 



z=0 



z=0 



are tenable in a < x < b, 16 K 1 , then Vx € [0, K^) 



Enirxo rnrx 1 u \ 1 viu; u \ . . 

a n cos cos = cosh hco-^— cos — — o(e ) 



n=0 



E fl ™ 

n=0 



d 



dx 



tx d 



nirxo . nirx ( d , 

cos sin = cosh xo — sin 

c c \ ax 



c <9z 

7TX d 

c dz 



S{e z 



z=0 



2 = 



(47) 



are tenable in a + xo < x < b — x$. Accordingly, for any definite value x € (a, 6), xo in (|47[) takes 
values in the following interval: 



< x < 



x — a, a < x < (a + b)/2, 
b- x, b> x > (a + 6)/2. 
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Where the left side of g7]) is: 



nirxo nnx ( d \ nitx 

cos cos = cosh x — — cos , 

c c V ox J c 



nnxo nnx ( d \ nnx 

cos sin = cosh xott~ sin . 

c c V ox J c 



Therefore, wxh(hd x ) and coah(hd x ) are the continuous operators on L 2 {[— c, c]). 

2.4 Abstract operators and Fourier analysis 

The third way to extend the domains and ranges of abstract operators is using Fourier analysis. 
Theorem 5. Let t G M\, V/(x) G ^(R 1 ), if F+(s), F_(s) are respectively 



where Jzf is the Laplace transform operator, then we have 

d 



f(x) = cos ( x— ) F+(s) 



s=0 



s=0 



Or 



cos (x-) F + (s) ^ 
Proof. By (jTTJ) we have 



f(x)+f(-x) 



d \ 1 



f(x) cos(tx)dx = ^f- 1 F+(s). 

X) 

/(a;) sin(te)dx = 5£- x F^ (s). 



cos(i£) = cos \i— )-, JSf sin(tf) = sin -, i G Mi, SR(s) > 

\ us J s \ os J s T 

So by (gU and (09}, V/(i) G ^(M 1 ), 3/ a (a:) G C 00 ^ 1 ) can be expressed as 



(48) 

(49) 

(50) 
(51) 



co> I x-^- I F + (s) + sin [ — x-^- I F_(s) 



cos a; 



as 



2 /"°° 

7T 



<9s 

/(Ojsfcos(*e)de 







+ sin ^ 





2 /"OO 
7T 



i / a 

— cos a;— 
7T \ as 

^ /*°° 

7T 



sin — x 



9s 



/(0 



J-oo V OS/ S 7T 

cos (x£) cos U£\ + sin sin 



/(OJSfsin(^)de 



s 



7T J-oo V OS/ S 7T J-, 

Therefore, we get 

lim / s (z) = lim - ^ f(0 

s^0+ s->0+ 7T J_ 00 



(x - e) 2 + s 2 

s 



/(0 



(x - if + s 



: di, K(s) > 0. 



fd)S(x - Odi = f(x). 
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In the same manner we find 

Theorem 6. Let f(x) € L 2 ([-l, I]), f(x + 21) = fix). If S+(t) and SL(t) arc respectively 



S+(t) = 777 /(0 



2/ J_i ^ J 1 -2tcos«//) +< 2 



Sl(*) 



' ./ ^ fsin«/Z) 

, - 2tcos«/0 +t 2 



then 3/ z (x) € C 00 ^ 1 ) with the following form 







/,(*) = cos — — S+(e z ) + sin — — 5_(e z ), -co < z < 



Z dz 



making 



lim f z (x) - f(x), Vf(x) e C(-Z,Z) C £ 2 (H,Z]). 
lim 11/(0;) - /,(a:) || i2([ _ M]) = 0, Vf(x) e L 2 ([-l, I}). 



Proof. By the algorithms ([21]) we have 

'< 9 \ 1 + e 2 



cos 



I dz 1 



1 -e 
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So / z (x) can be expressed as 

/ 7TX d 

/,(*) = cos -- 



1 - 2e z cos(tt£/Z) + e 2z ' 

e z sin(7r£/Z) 
/ 9z7 1 - e z ~ l-2e z cos(tt£/Z) + e 2z ' 

?r£ 9 \ 1 + > ■ 



f($ COB 



-/ 



Sill 



1 

2Z 



7TX 8 

~T~d~z 

f(t) cos 



I dz I 



tt£ d 



:dt 



-! 



/(Osin 



7TX 5 
7TX 8 

Tdz 



I dzj l-e z 
2e 



dt 



1 



Tttd_ 

~T~d~z 
TdzJ 1 



1 



:dt 



f(0 



cos 



TTX 8 

~T~d~z 

I 



COS 



777 /(£K+7 /(Ocos 



n(x-g) d 



sm 



2/ 

7rir d 
~T~dz 



fiOdt 



-/ 
sin 



<_9 
Z <9z 



Z 9^/1 



27 7 /(0 - 
7TT /' /(0 



?r(x - 9 \ 1 



9z / 1 — e 

■ e 



2/ vsy 1 - 2e z cos(7r(a; - £)/Z) + e 2 



-oo < z < 0. 



(52) 
(53) 

(54) 



7# 



1G 



Therefore, if f{x) G C(-l, I) C L 2 ([-l, /]), we get 



1 f l 

hm f z (x) = lim — / /(£)-; — 



1-e 2 



/(C) Hm 



1 



2e z cos(tt(x - £)/0 + e2z 
l-e 2z 



_ t z^o- 21 1 - 2e z cos(tt(x - f)/0 + e 22 
As we all know, V/(x) € L 2 (H, J]), /(x + 21) = f(x) we have 



f(Q6(x-t)dt = f(x). 



,. ^ / knx , . fc7ra;\ 
/ (x) = a + I Ofe cos — 1- 6fe sin — — I 



fc=i 



Where Vfc G Nj 



li f /(Ode, a* = j /' /(e) cos 6 fc = 1 f l m) sm 



Accordingly f z (x) G C°° (M ) can be expressed as 

f z (x) = a + 22 ( a fe efcz cos —j 1" bk£ kz sin — — ) , -oo < z < 0. 



fe=i 



So we have 



- / \f(x) - /.(x)! 2 ^ = £[(a 2 + 6 2 )(1 " e fez ) 2 ], -co < z < 0. 

fc=i 

lim ||/(a:) -/,(a;)|| il(H ,q) = 0, Vf(x)eL 2 ([-l 7 l}), f(x + 2l) = f(x). 
If /(cc) is an even function, then we have the integral equation: 



1 



1-t 2 



I J 1 - 2tcos(irx/l) + t 2 
Its solution f(x) G L 2 ([— I, I]), f(x + 21) — f(x), and we have 

f{x)=m J™d-\ S+ (e*) 



dx = S+{t), VS+(t) G C°°(O ). 



z=0 



If f(x) is an odd function, then we have the integral equation 

t sm{TTx/l) 



2 

T/o 



1 - 2tcos(nx/l) +t 2 



dx = S-(t), yS-(t) G C°°(fio)- 



Its solution /(a;) £ L 2 ([-(,l]), f(x + 2l) — f(x), and we have 



/(a;) = sin 



7TX d 

Tdz 



SAe z ) 



(55) 

(56) 
(57) 

(58) 
(59) 
(60) 
(61) 



For instance, in the integral equation (|58p. if S+(t) = arctani, < t < 1, then its solution is a 
square wave, namely 



/(*) 



+tt/4, 2Kl-l/2<x <2Kl + l/2 
-tt/4, 2KI + 1/2 < x < 2KI + 31/2. 
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Where K = 0, ±1, ±2, • • •. 

Let p z (x) £ Cg°(R ra ), svppp z (x) C {a; £ R n | |x| < 1} be a polished kernel, taken as p z {x) — 
p Zl (xi)p Z2 (xs) ■ ■ ■ p Zi (xi) ■■■ p z „ (x n ) and 

pM = Wa-2e«™{«x-/l i )+e*>'> -~<*.<0,N<« i (i = l 1 2 | ...,n). (62) 
Then we easily prove the following results: 

hm g{d x ) [ f{t)p t (x-Z)dt^g(d x )f(x), V/(z) G L a ([-U]), .9(0 G C°°(R n ). (63) 

Where /(x + 21) = f(x), I £ W'\ [-1, 1} = {x £ R n | ~l<x<l}. 

If g(d x ) — d x , Va G N", then (|63l) gives a definition of generalized derivative on L 2 ([—l,l]). 



2.5 Integral representation of abstract operators 



Clearly, such algorithms as (j20|) - (f25|) do not exist for most abstract operators. Therefore, 
for more complex abstract operators, establishing the integral expression is the fourth important 
way to extend the domain of abstract operators. In reference [2j Guangqing Bi has obtained the 
following results: 

Theorem BI1. Let P(d x ) be an m-order partial differential operator of any kind and there 
exist ai, a2, • • • , cik of real and partial differential operators A\, A2, • • • , A^ of the order less than 

[(m + l)/2] such that P(d x ) = a x A\ + a 2 A\ H V a k A\. If k = 2v + 3, v = 0, 1, 2, • • •, then 

Vf(x) £ C mv (R n ) we have 

_ t r\ dt ...r tdt (^,)r 




fc-2 
f 2i+l 



E^TlM^^'/C*)- (64) 



(2 

i=0 v 

Where a is a real numbers, a; G R™, t G R 1 . 7? G Rfe is the integral variable and 

rj\ = at cos 0i , 

r/2 = at sin 6\ cos #2, 



f? p = at sin sin 02 • • • sin0 p _i cos9 p , 
rjp+i — at sin 0i sin 02 • • • sin p cos </>, 
rjp+2 — rjk — at sin 0i sin 02 • • • sin p sin (j); 
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da k = sin k 2 6\ sin fc 3 9 2 ■ ■ ■ smd k -2d8 1 d92 ■ ■ ■ dOk-idtj). 
When Ai, ■ ■■ , Af. in the right-hand of are one order partial differential operators, then the 

1/2 A , , 1/2 A 

abstract operator e ??ia i lH A k j g Qne Q £ ^ e following five simplest operators: 

exp(hd x ), sin(hd x ), cos(hd x ), smh(hd x ) and cosh(hd x ). 
Proof. In (JHll), let f(x) = e^ x , x G R n , £ £ R n , and the symbol of the partial differential 

1 /2 

operator Aj,j = 1, 2, • • ■ , fc is denoted by Xj(£)> ft = a - Xj(0- By Definition 4, the formula 
degenerates to its characteristic equation: 

(a «p( {) v=) _ tdt yp(i))- r r.„ r e „^ Adai 



f fe-2 

^ t 2i+1 k - s 

i=0 v ' 

According to the analytic continuous fundamental theorem, we only need to prove (|65p . Solving 
the integral on a hypersphere on the right side of (|65[) . we have 

sinh(atP(0 1/2 ) _ /■* f'tjt'Sr (a 2 P(0)" +j t 2j ^ t 2l+1 (a 2 P(0) J 



^K) 1/2 7o 7o £j(2j)!!(2j + 2i/+l)l! + ^ (2i+l)l ' 

Then it is proved by the termwise integration of the infinite series on the right side of the equation. 
Similarly, we have (See [1]) 

Theorem BI1'. Let P(d x ) be an m-order partial differential operator of any kind and there 
exist ai, a%, ■ ■ ■ , of real and partial differential operators Ax, A 2 , ■ ■ ■ , Af. of the order less than 
[(m + l)/2] such that P(9 X ) = aiA 2 + a 2 v4| H h a fc A^. If k = 2v + 3, v = 0, 1, 2,- ••, then 

V/(z) g C°°(C n ) 



* f tdt-- - [ tdt[-P{d x )} u 
Jo Jo 



v+l r n/2 ptt/2 

V2 A \ v„ 1/2 



[~J I "'J cos{r] 1 a 1 1 / ^A 1 )---cos(r] k a L k / ^A k )f(x)da k 
fe-i 

E^TilTl-mF/W. (66) 



(2i + l)l 

Where x G K", t 6 I 1 . 77 G Kfc is the integral variable and 

rji = t cos 9 1 , 

rj2 — t sin #1 cos #2 , 
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r]p = t sin 0i sin 02 ■ ■ ■ sin0 p _i cos0 p , 
r) p+ i — t sin 6\ sin 82 ■ ■ ■ sin p cos 0, 
7?p+2 = Vk — t sin 9\ sin 02 ■ ■ ■ sin 6 p sin 0; 

da k = sin fc_2 #i sin fc ~ 3 9 2 ■ ■ ■ sm9 k ~2d9id9 2 ■ ■ ■ d9 k -2d(j). 

If A\, ■ ■ ■ , Af. in (|64p and (|66l) are partial differential operators of the order great than 1, then 
the order can be lowered by taking the following Theorem: 

Theorem BI2. Let P(d x ), x G M. n be a partial differential operator of any order and f{x) G 
C(M") which make the integral in the follow formula meaningful, then 

e xp(c%) f{x) = _J_ / e -t 2 /^ 1/2 (Hd*)^ f{x)dC: VA G C 1 . (67) 

2V7T , 



Continuing the processing one can arrive at order one. 

Example 6. By using the Theorem BI2 and Corollary 2 we have 



cxp [ —a 2 t-^-^ I q(x) = \= I e ^ l A cos ( a 
y \ dx] J yv ; 20F 



V~t(-^g(x)d(, t€R\. 



Where ^,o(C) = aA + yO^C, VA,a e C 1 . 

With the expansion of definitions of abstract operators, it is necessary to add new algorithms 
continuously. 

3 Abstract operators and linear higher-order partial differ- 
ential equations 

3.1 Applications of Laplace transform for n-dimensional space 

Solving the ordinary or partial differential equations, is constructing the algorithms of the 
inverse operators of ordinary or partial differential operators. In terms of abstract operators, 
solving the initial value problem of partial differential equations for n-dimensional space is similar 
to solving the ordinary differential equations with respect to the variable t, thus we can introduce 
the Laplace transform to further simplify the solving process. 

Corollary 7. Let g(J*-tdt) be the abstract operators taking J^-tdt as the operator element, 
and g(— ijj) be the abstract operator taking — ^-3- as the operator element, and the symbols of 
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these two abstract operators are the same, then for any Laplace transformable function f(t) we 
have t 

g[J-lM^f(t)=S>- 1 g(-~^Sff{t), teW^. (68) 

Where Jzf and _Sf _1 are the Laplace transform operator and its inverse operator respectively, and 
J?f(t)=F(s). 

Proof. According to the properties of Laplace transform, we generally have 
J -tdt) f(t)=J?~ 1 {-- — J F(s), VmeN. 

According to the abstract operator fundamental theorem, we get the Corollary 7. 
Similarly we have 

g(J*-dt^f(t)=J?- 1 g(^)j?f(t), t€M\. (69) 

Therefore, using properties of the Laplace transform and the abstract operator fundamental theo- 
rem, we respectively obtain the abstract operators taking f*- tdt and J*- dt as the operator elements. 
This type of formulas can also be applied in solving certain integral equations. 

Theorem 7. Let m € N+, P{d x ) be arbitrary order partial differential equations for n- 
dimcnsional space. Then we have 



d 2 

0^-P(d x )) u = f{x,t), iel",fet|,VmeN + , 

&>u 

~W u-, 



<Pj(x), j = 0,1,2,..., 2m- 1. 



(70) 



U{X ' t) = JoJ (2m -2)1! (2m -4)1! P^ 2 f^)r dr dr 

m ~ 1 / m \ 2m-l-2fc Q2m-l-2k-j rt / f 2 _ r 2-jm-2 T 

+ £ ( k \k) P{ - dx)k ^ dt 2 ™- 1 -™-' J (2m -2)!! (2m -4)!! 
sinhfrP^) 1 / 2 ) , N 

Proof. The Laplace transform of the partial differential equation, with respect to t and con- 
sidering the initial condition, is 

m , / 2ro-l-2fc \ 

Where U(x,s) = 3?u(x,t), F(x,s) = S?f(x,t). Let G m (d x ,t) = jSf" 1 ^/^ 2 - P{d x )) m ], solving 
U(x, s) and its inverse Laplace transform is 

u(x,t) = ^U(x, S )=J?- ^ s2 _^ ))m F( X , S ) 
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m ~ 1 /™\ 2m-l-2ft „2m-l-2fc-j 

G m (d x ,t)*f(x,t) 



2m-l-2fc 

m 11 



, Q2m— 1 — 2fe— j 

Q^2m-l-2k-j Gm{dx,t)<Pj{x). 

Now let us solve G m (d x ,t). In Corollary 7, taking = £ m_1 , £ G Mi as the symbol of the 
abstract operator, and let f(t) = sin&i, t G K 1 , we have 



Let b = iP(0 1/2 , £ S R„, we have 



1 _ 1 / /■* . ,A m_1 sinh(tP(QV 2 ) 



: = 77. ~T~77 [ / -tdt 



(s 2 -P(0) m {2m-2)\\\J J P(0 1/2 

Taking this one as the characteristic equation, according to the analytic continuous fundamental 
theorem, we have 

, s i 1 Iff* A™" 1 swh(tP(d x ) x / 2 ) , 



■xdx) f(x)= xdx--- f(x)xdx= \ n s /(flfrg- (73) 



By using (8) in reference [3], we can easily derive the following integral formula 

r (z 2 - C 2 ) m " 
/„ (2m -2)!! 

Applying ([73)) to (|72|) . we have the expression of abstract operator G m (d x ,t): 

r* ( t 2 _ T /2)m-2 smh ^Jp^l/2\ 



G m (d x ,t)g(x) = J q (2m : 2)!! ( 2m „ 4 )!! pfcp? ' 9 {x)r'dr'. (74) 

Thus Theorem 7 is proved. 

In reference [5] the Guangqing Bi has obtained the following results: 

Theorem BI3. Let ai, a®, . . . , a m be arbitrary real or complex numbers different from each 
other, P{d x ) be a partial differential operator of any order, then we have 



(75) 

= 0, j = 0,1,2,...,™-!. 



m „ 

TT(-- 






~dV~ 


t=0 




r 


lo J 






, m — 1 



/, /\m — 2 '"• /I 

««'■'»-/./. ' ( ;V 2 j| r ng ,^^,^'/(x.r)^. (76) 
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Theorem BI4. Let a±, ■ ■ ■ , a m be arbitrary real or complex numbers different from each 
other, P{d x ) be a partial differential operator of any order, then we have 



d 



(77) 



= 0, 



j = 0,l,2,...,2ro-l. 



u(x, t) = 



t rt-T n _ -/\2m-3 



{t-r-r'Y 



E 



,2ro-2 



smhir'ajPidx) 1 / 2 ) 



(2m 3)! ^ n™- K 2 - ) a^)^ 



f(x,T)dT'dT. (78) 



On this basis, by using the abstract operators and Laplace transform we have obtained the 
following results: 

Theorem 8. Let ai,a 2 , ■ ■ ■ ,a m be arbitrary real or complex roots different from each other 
for bo + bix + &2X 2 + • • • + b m X m = 0, and P(d X) d t ) be a partial differential operator defined by 



fc=0 



"' 0j*> xeR n 7 tem. 1 + ,VmeN + . 



Where P{d x ) is a partial differential operator of any order. Then we have 



P(d x ,d t )u = f(x,t), x e K", t <G R]_, Vm e N+, 
Q r 11 

= ip r (x), V = 0,1,2,..., 772 — 1. 

t=0 



dt r 



(79) 



u(x, t) 



(t-r-r') 



f\m-2 m 



o j (777-2)! 4^n^(oj-oi) 



E 



SajPiO.)^^ dT > dj . 



' 1 Qk — l-r j-t _ T \m-2 



k=l r=0 



(m-2)! 



E 



— 1 



</? r (a:) dr. 



(80) 



Theorem 9. Let a\, a 2 , ■ ■ ■ , a m be arbitrary real or complex roots different from each other, 
satisfy X^fcLo ^2kX 2k — ll"=i(x 2 ~ a l)> an d P(d x ,d t ) be a partial differential operators defined by 

m f)2k 

m^d^^mr^, xgm", teR\, V777 GN+. 



k=0 



Where P(d x ) be a partial differential operator of any order, then we have 



P{d x ,d t )u = f(x,t), iel",ie R\, Vm e N+, 

= <p r (x), r = 0,1,2, ... ,2t77 — 1. 

t=0 



9r 



(81) 
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(,J kh (2m -3)1 f^n™-K 2 -af) 0j P(fi!.)V» 

(«?-«?) ^m) 1 / 2 M) (8) 

Let us prove these two theorems. According to the Theorem BI3 and Theorem BI4, we just 
need to prove the following Corollary of Theorem 8 and Theorem 9: 

Corollary 8. Let a±,a2, ■ ■ ■ ,a m be arbitrary real or complex roots different from each other 
for b + b\X + &2X 2 + ■ ■ • + b m x m = 0, and P(d x ,d t ) be a partial differential operators defined by 

P(d x , d t ) = b kP(d x ) m - k Qj:, xeR n ,te R\, Vm e N+. 

Where P(d x ) be a partial differential operator of any order, then we have 



P(d x , d t )u = 0, x e M n , t € M.I, Vm e N+, 

= <p r (x), r = 0,1,2, ...,m- 1. 

4=0 



9r 



(83) 



**> t] = £ fc^r-* Etp^/ T^yr E fc^y ^ dT - (84) 



fc=l r=0 



Proof. Considering initial conditions, the Laplace transform of the Eq ((53)) with respect to f 

is 

m / fe-1 \ 

E W^)" 1 "* *) - E S k - X - r ?r{x) = 0. 

fc=0 V r=0 / 

Where C/(x, s) = J?u(x,t), considering IIi=i( s — a iP(9 x )) = Sfclo bkS k P(d x ) m ~ k we have 

m m ft — 1 

H(s - ai P(d x ))U(x, *) - E W^)"* - * E s fe_1 ~V(*) = 0. 

i=l fc=l r=0 

We need to introduce an abstract operators G m (d x ,t) 7 defined as 

G m (d x , t) = J£ 1 - 



n£i(»-aiP(Sx))' 
By solving U(x, s), we have its inverse Laplace transform: 

m fe — 1 fc _ i 

k— 1 r— IV J V a:// 

7Tt fc— 1 r\fc_\_ r 

= Y, b ^ d -) m ~ k Y,w^ Gm{d ^ t)lfr{ - x) - (85) 

fc=l r=0 
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Now let us solve G m (d x ,t). Considering initial conditions, the Laplace transform of the Eq 
f5|) with respect to t is 

Y[(s- ai P(d x ))U(x, s ) =F(x,s), 

i=l 

where F(x,s) — J£f(x,t). By solving U(x,s) and using the convolution theorem, we have its 
inverse Laplace transform: 

u(x, t) = %- x U{x, s) = if" 1 1 F(x, s) = G m (d x ,t) * f(x, t). 

n l= i( s ~«^(^)) 

By comparing (|76|) with u(x,t) = G m (d x ,t) * f(x,t), we have the expression of the abstract 
operators G m (d x ,t): 

G m {d Xl t)= { ) T} J2fr^j -^ P ^dr. (86) 

Jo (m -2)! ^ n*=H ffl j 

Applying ([8"6")) to (|8"5)l . thus the Corollary 8 is proved. 

Corollary 9. Let a\, a-x, . . . , a m be arbitrary real or complex roots different from each other, 
which satisfy X^fcLo ^2fcX 2fc = YITLi (x 2 — a f X an d P(d x , d t ) be a partial differential operator defined 

by 

p(d x ,d t ) = J2 b ^ p ^ m ^ k -Q^k^ X G R "' 1 £ R +' Vm g n +- 

Where P{d x ) is a partial differential operator of any order, then we have 



P(d x , d t )u = 0, x G K", i G R^, Vm e N+, 

d r u 

— (pr(x), r = 0, 1, 2, . . . , 2m — 1. 



-' 1 Q2k-\-r r-t _ T j2m-3 



(87) 



^ 5t 2fe - 1 -'' 7 (2m- 3)! 

k=l r=0 Jy> \ ) 



2m-2 

X 

3=i 



\- «T : sinh(r aj P(^) 1 / 2 ) 

1 i^J ■> 



Proof. Considering initial conditions, the Laplace transform of the Eq (|8"T)) with respect to t 

is 

m / 2k-l \ 

J2 b ^P(d x ) m - k is 2k U(x,s)- £ ^- 1 -> r ( 2 ;) =0. 

fe=0 V r=0 / 

Where C/(x, s) = Sfu(x, t), considering J]™ i(s 2 - afP{d x )) = J2T=o hkS 2k P{d x ) m - k we have 

m m 2/e— 1 

Il(s 2 - a 2 P(^))(7(x, s)-J2 b2kP(d x ) rn - k £ s^^VM = 0. 

i=l k=l r=0 

We need to introduce an abstract operators G m (d x ,t), defined as 

f 



G m (d x , f) — 



Il£i(« a -a?JW)" 
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By solving U(x, s), we have its inverse Laplace transform: 



2fc-l 



c(x,t) = X- l u(x,s) = J2b2 k P(d x ) m - k E-^ fF (2- *Pm w V'to 



k=l 



„ nr^ 2 -^*)) 

m 2 A:— 1 o2fc— 1 — r 

E^m)™^ E ft2fc-l-r G m {d x ,t) Vr {x). (89) 



Now let us solve the G m (d x , t). Considering initial conditions, the Laplace transform of the Eq 
(1771) with respect to t is 

m 

l[{a 2 -a$P{d x ))U(x,a) =F(x,s), 
»=i 

where F(x,s) = Jzff(x,t). By solving C/(x, s) and using the convolution theorem, we have its 
inverse Laplace transform: 

u{x,t)=J?- 1 U{x, S )=J?- X Tim f2 1 T = n r^r F(x, S )=G m (d x ,t)*f(x,t). 

By comparing (|75|) with u(x,i) — G m (d x ,t) * f(x,t), we have the expression of the abstract 
operators G m (d x ,t): 

G (d t)= [ Ht-r^-' f a*- 2 sinh(, flj P(^) 

Applying (|9"0")l to (f8"9")l . thus the Corollary 9 is proved. 

3.2 Analytic solution of Cauchy problem 

Theorem 10. Let A„ = + + • • • + g§-y be an n-dimensional Laplacian. If n — 2 = 
2i/ + 1, i/ eN, then V/(af,r) £ C 2,y (E"), and V^(a;) G C 2m+2 ' y - 2 (M"), we have 

— -a 2 A„) u = f(x,t), x G M™, t G K^, Vm G N+, 

d3 



/ (91) 
^■(x), i = 0,1,2,..., 2m- 1. 



u{x,t) = 



, ra-2 



7 (2m -2)!! (2m -4)!! 

r ' L2A M/ 



'2 



32)!! ? I (2 Ji + 2r-l)!l(2r)I! ^ r ^ r ^ 



(2m 



r=0 ' 



2m-l-2fc Q2m-l-2fe-j /-t / + 2 _2\m-2^2 



k=0 



Qfim-i-sk-j J Q (2m - 2)!! (2m - 4)!! 
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2m-l-2fc 



" 1 'm-l + A t2fc+2i+j 



X 

i=0 



Where n- 2 = 2i/ + 1, S' n = 2(27r) l/+1 (ar')"" 1 , 5„ = 2(2vr) l '+ 1 (ar)™- 1 , and f € M„ is the integral 

variable. The integral is on the hypersphere (g[ — xi) 2 + (£' 2 — x 2 ) 2 H h — x„) 2 = (ar') 2 , and 

dS^ is its surface element. £ € R„ is the integral variable on the hypersphere (£i — xi) 2 + (£2 — 
X2) 2 + • • • + (£ n -i n ) 2 = (ar) 2 , and dSVi is its surface element. 
Proof. According to (j64]), we can easily derive: 

sinh(atA„ 1 / 2 ) /s f* f* (a 2 A n Y f 

'-f{x) = t tdt--- 1 " j / f(OdS n tdt 



GlA^ 1 / 2 Jn Jt) S n J g n 



^ t 2i+1 

+ E nr(« 2 A„)7W, V/(z) e C 2 "(R"). (93) 

Where A„ is an n-dimensional Laplacian, n — 2 = 2iv + 1, S n = 2(27r) ly+1 (ai)™~ 1 . £ e R„ is the 
integral variable on the hypersphere (£1 — x\) 2 + (£2 — x 2 ) 2 + ■ ■ ■ + (£„ — x„) 2 = (at) 2 , and cZS n is 
its surface element. 

In Theorem 7, let P(d x ) = a 2 A n , then Theorem 10 is proved by the substitution of 
Theorem 11. Let a\, a 2 , ■ . . , a m be arbitrary real roots different from each other, which satisfy 
Y^k=o^kX 2k — Il"=i(x 2 — a f)i an d P{d x ,d t ) be a partial differential operator defined by 

m a2k 

P(d x ,d t ) = & 2feA™- fc ^, xeR n ,te R\, Vm e N+. 

fc=0 

Where A„ is an n-dimensional Laplacian. If n — 2 = 1v + 1, ^ = 0,1,2, then \/ip r (x) £ 
C 2m+2l/ - 2 (R"), and V/(x,t) € C 2l/ (R"), we have 

P(d x ,d t )u = f(x,t), x e R n , t e iy, Vm e N+, 
<9 r tt 



9#" 



(f r (x), r = 0, 1, 2, . . . , 2m — 1. ^ 



t=o 



u(a;, f) = 



/0 (2m -3)! ^nf-K 2 - a ? 



' f r'dr' ■■■[' / /(£', r) dS' n , r'dr' + g 



,-i a2 ; r ,2i+i 

- — Aj,/(x,r) 



(21 + 1)! 



2 — 1 Q 2k - l -' r '"' 



fe=l 



r=0 



^2fe-l- 



r/r 



(t-r) 



2m-3 



E 



,2m— 2 



(2m -3)! ^n^(o»-a?) 
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r T (a 2 A Y r V ~ X a 2l r 2l+1 

Tdr . . . I ±J-2L- / <p r (£) dS n<j TdT + V -1-—A l n ip r (x) 



\ — : — / 



(95) 

Where S' n ] = 2(27r) I/ + 1 (a 3 T')"~ 1 , S nJ = 2(2 7 r) I/ + 1 (a i r)"- 1 , and £' € R„ is the integral variable. 
The integral is on the hypersphere (£j — xi) 2 + (£ 2 — a^) 2 + ■ ■ ■ + (£ n — x n) 2 — ( a j T ') 2 , and dS' n j 
is its surface element. £ £ R„ is the integral variable on the hypersphere (£1 — xi) 2 + (£2 — a^) 2 + 
■ ' ' + (£n — a;n) 2 = ( a j T ) 2 j and dS'nj' is its surface element. 

Proof. In Theorem 9, let P(d x ) = A„, then Theorem 11 is proved by the substitution of (|9"5]l . 

Similarly, we can easily obtain explicit solutions of the Cauchy problem of more complex partial 
differential equations. Therefore, using the method of abstract operators, we have established a 
general theory of initial value problems for linear higher-order partial differential equations. 

3.3 Initial-boundary value problem and Hilbert space 

Clearly, similar to smh(hd x ) and cosh(hd x ), the abstract operators 

sinh(a<A„ 1 /2) q sinh A n x / 2 ) 

aA Tt V2 and COs r A " ) = di Q A„V2 

on the left side of (|93p are bounded operators in a Hilbert space, and can act on the whole Hilbert 
space. Thus we can attach proper boundary conditions to the initial value problem introduced 
by Theorem 7, which makes the given functions f(x,t), <fj(x) in Eq (|91[) and (|94[) become func- 
tions with boundary conditions. Then solving this initial-boundary value problem comes down to 
expressing f(x,t), <Pj{x) in a Hilbert space H within the given domain Q. 

For the initial-boundary value problem, the operator P{d x ) must have the characteristic func- 
tion related to boundary conditions, in order to expand the known function f(x, r), tp r (x) in (|71[) , 
(1501 ) and (|82p by using the characteristic function of P{d x ). P(d x ) in Theorem 7, Theorem 8 and 
Theorem 9 can be variable-coefficient partial differential operators. For instance, if P(d x ) is a 
self-adjoint operator defined on a Hilbert space H, then the abstract operators 

sinhftP^V 2 ) , _ 9 sinhftP^V 2 ) 

P(0 a )Va and C0S *H^W J ^ P(0 a )Va 

are the continuous operators on the Hilbert space. In this case, we can attach proper bound- 
ary conditions to the initial value problems in ([70)1 . (|T9"|) and (|8ip . Therefore, the given function 
f(x,t),tp r (x) becomes a function with boundary conditions, and can be expressed in the Hilbert 
space H within the given domain CI. In order to solve the corresponding initial-boundary value 
problem, we need to solve the characteristic value problem of P(d x ) under given boundary condi- 
tions to determine a set of orthogonal functions, which generates a linear manifold of the Hilbert 
space, thus f(x,t),(p r (x) can be expressed in the Hilbert space. 

In (|TTj) . (|5U1) and (|B"2"|) . if f(x, r), ip r (x) € L 2 (Q), and P(d x ) is a second-order linear self-adjoint 
elliptic operator, namely 

n 8 / (97/ \ 

P(3r)ii = V t— Oij(x)^- + c{x)u, xeilc R'\ (96) 



2cS 



then boundary conditions can be added for the definite solution problems (|70p . (|79p and (fSTj) : 
Bu\gn = Orepresenting = or 



a,ij(x)-j^- cos(a,Xi) + 6(a;)t 



= 0. 



<■«> 



Where a is the unit outward normal of dfl. Thus this kind of initial boundary value problems 
boils down to solving the characteristic value problem of the first boundary value problem of 
second-order linear self-adjoint elliptic operator: 



E 

i,j=l 



d_ 

dxi 



dij(x)— — ) + c(x)u = — Aii, ie!!c 

OXi 



(97) 



Bu 



an 







Reference [7] quotes the following theorem: 

Let fl C K n be a bounded open domain, and dQ be smooth. Let ay — dji, and there exists 
6 > such that 

n 

0|£l 2 > xeCi. 

And let ay G C^H), c(x) e C(H), i>(a:) G C(<9f2), then ^7} has the following countable charac- 
teristic values: 

< Ai < A2 < ■ ■ • < X v < ■ ■ ■ , lim A„ = 00 

v— >oo 

(If (ay) = J is a unit matrix, then Ai = when b(x) = c(x) = 0. When 6(2;) > 0, c(x) > and one 
of them does not identically equal to zero, Ai > 0) and the corresponding characteristic functions 
ei(x),e 2 (x), ■ ■ ■ ,e v (x), • • • , satisfy 



E 



d_ 

dxj 



d& ■ 

a ij( x )~Q^; ) + c ( x ) e i = -AjSj, (e,, e,) = (5j 



(98) 



and {ej(x)}jL 1 is complete in L 2 (f2), that is for an arbitrary f(x) G L 2 (f2), there exists <y such 
that 

m 

iim 11/- E c J e jii L2 ( o ) = °- 

3=1 

Therefore, for ([7T]). ([50)1 and ([521), by f(x,r),(p r (x) G L 2 (fJ) we also have 

m 

lim ||/(2, t) - Vc 3 (r) ej (a;)|| L 2 (f2) =0. 
3=1 



lim ||y r (ai) - E c J' e J'( x )lli 2 (") = °- 

m— )-oo * — * 
3=1 

Clearly, based on ((96|) and (l98l) . we have the extension of abstract operators in Hilbert space: 



_„ P /a \ , . >, , . sinh(ra ? -P(9 a; ) 1 / 2 ) , . sin(a 7 --v/A7 r) . . . . 



Q>i V Ai 
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Where x £ Q, C E n , j = 1, 2, • • • , m. Then we have 



lim \\e^ p ^f(x,r)~J2^ e ^ T ' a3X ' e ^\\^ 



lim 

m— ^oo 



lim 

m— >oo 



hm ||e ro * F ^Vr(a:)-y](Vr > e i )e- ra * A 'e i (a!)|| za( n)=0. 



sinh(r'a j P(^) 1/2 ) -A sin(%'VAlr') 



j'=i 



sinh(ra i F(5 x ) 1/2 ) sin^VAiV) / N 

Oj .P0,)1/2 yr(x)-Xj( yr , ej ^= e,(») 



L 2 (H) 



L 2 (£l) 



Thus when (f7U|) . (|79p and (j8"Tj) satisfy the boundary condition Bu\qq_ — 0, if the orthogonal 
function system {ej(x)}j2. 1 is known, we have the solution of the corresponding initial boundary 
value problem. 

In conclusion, based on the method of abstract operators and comprehensively applying the 
Laplace transform and Hilbert space theory, we establish a general theory of initial value and initial 
boundary value problems of linear partial differential equations, and illustrate broad application 
prospects of the abstract operator theory, which will definitely have a wide and profound influence 
on various scientific fields. 
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